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SECTION  I 


INTRODUCTION 

A.  PROBLEM,  OBJECTIVES,  AND  PURPOSE 

As  part  of  a current  crash  avoidance  research  initiative,  NHTSA  has 
been  looking  into  the  relative  occurance  of  various  types  of  accidents. 
When  car  size  is  taken  into  account,  it  has  been  found  that  the  likeli- 
hood of  roll  over  increases  with  decreasing  car  size.  In  an  effort  to 
uncover  factors  that  might  contribute  to  this  accident  in  lvement  mode, 
NHTSA  has  been  working  with  a relatively  simple  model  for  describing 
tripped  roll  overs  (Ref.  1).  This  model  accounts  for  a vehicle  sliding 
sideways  into  a curb  and  can  predict  roll  over  potential  based  on  ini- 
tial side  velocity,  and  certain  details  associated  with  the  location  of 
the  vehicle  center  of  gravity  and  roll  characteristics  of  the  vehicle’s 
sprung  mass.  This  modeling  approach  appears  to  have  some  potential  for 
expanding  our  knowledge  of  the  factors  contributing  to  vehicle  roll  over 
involvement  in  accidents,  but  the  generality  of  the  model  must  be 
expanded  in  order  to  fully  explore  this  method. 

The  Ref.  1 vehicle/curb  encounter  model  is  fairly  restricted  to 
vehicle  velocities  normal  to  the  curb,  with  both  wheels  encountering  the 
curb  simultaneously.  The  objective  of  the  work  reported  herein  was  to 
develop  an  expanded  computer  model  to  include  several  arbitrary  initial 
conditions:  1)  curb  encounter  angle;  2)  vehicle  lateral  velocity; 
3)  vehicle  yawing  rate.  The  time/force  relationship  for  the  wheel/curb 
impact  was  also  to  be  modelled.  The  Ref.  1 model  assumes  complete  con- 
servation of  the  impulse  momentum  associated  with  wheel/curb  side 
impact.  Some  consideration  was  also  to  be  given  to  energy  dissipation 
and  storage  due  to  wheel/suspension  deformation  and  elastic  deflection. 
Elaborate  geometric/structural  models  were  not  desired  here,  but  the 
model  should  account  for  simple  effects  such  as  deformation  and  com- 
pliance, and  damping  associated  with  the  combined  elastic  deformation  of 
the  tires,  wheels,  and  suspension. 
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The  purpose  of  the  work  requested  under  this  task  order  is  to  pro- 
vide a vehicle  tripped  roll  over  model  that  NHTSA  can  use  to  further 
explore  factors  associated  with  this  accident  mode.  This  working  paper 
documents  the  implementation  of  the  model  in  a computer  program  which 
allows  initial  conditions  to  be  conveniently  manipulated,  and  provides 
for  convenient  graphical  output  of  results.  This  working  paper  docu- 
ments the  model  equations  (Appendix  A),  program  structure  and  prelimi- 
nary results  (Section  II),  and  includes  user  instructions  (Section  III). 

B.  BACKGROUND 


A fairly  simple  momentum  conservation  model  is  developed  in  Ref.  1 
to  determine  vehicle  roll  over  given  lateral  impact  with  a curb.  Given 
that  both  wheels  simultaneously  impact  the  curb,  the  critical  lateral 
velocity  for  overturning  is  given  by  the  equation 


where 


(1) 


g = gravity 

Mj,  = total  vehicle  mass; 


I = total  roll  moment 
of  inertia 


h„„  = center  of 

eg 


gravity  height; 


T = track  width 


It  is  clear  from  Eq.  1 that  the  critical  velocity  for  tip  over 
directly  decreases  with  the  ratio  of  center  of  gravity  height  to  one- 
half  of  the  track  width.  Therefore: 

hP£T 

Roll  Over  Tendency  = Y]2  (2) 

with  the  roll  over  tendency  increasing  with  decreasing  critical  veloc- 
ity. A survey  was  made  of  a range  of  domestic  and  foreign  cars,  and  the 
distribution  of  the  roll  over  tendency  parameter  is  plotted  in  Fig.  1. 
The  average  runs  about  0.75  for  domestic  cars,  while  the  average  is 
higher  (about  0.79)  for  foreign  cars  indicating  a potentially  higher 
roll  over  tendency  for  foreign  cars. 
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Figure  1.  Roll  Over  Tendency  Distributions  for  Domestic  and  Foreign  Cars 


The  conditions  leading  up  to  roll  over  accidents  are  not  very  well 
documented  in  accident  data  bases.  We  do  have  some  experimental  data 
from  a braking  experiment,  however  (Ref.  2),  where  the  vehicle  could  be 
configured  to  cause  the  rear  brakes  to  lockup  first  under  significant 
deceleration  conditions  (the  stopping  task  required  about  0.3  g deceler- 
ation). Rear  brake  lockup  can  lead  to  directional  instability  (Ref.  2) 
and  yaw  rate  was  recorded  for  all  test  runs.  The  distribution  of  yaw 
rate  was  analyzed  as  illustrated  in  Fig.  2.  Data  analysis  also  consid- 
ered the  covariation  of  heading  angle  and  speed  as  the  vehicle 
encroached  on  the  cone  delineated  lane  layed  out  for  the  braking  task. 
These  results  are  shown  in  Fig.  3.  The  data  in  Figs.  2 and  3 thus  gives 
some  indications  of  conditions  which  might  lead  to  tripped  roll  over. 

In  the  next  section,  we  will  describe  the  design  and  validation  of 
the  computer  program  which  was  prepared  to  analyze  the  the  roll  over 
tendencies  of  conditions  such  as  illustrated  in  Figs.  2 and  3. 
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Distribution  of  Absolute  Peak  Yaw  Rate  During  Braking  Tests 
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SECTION  II 


PROGRAM  DESCRIPTION 

A.  DEGREES-OF-FREEDOM  AND  DYNAMIC  MODES 

The  computer  programs  developed  during  this  project  simulate  a car 
sliding  into  a curb  using  a multi  degree-of -f reedom  model  and  various 
dynamic  modes.  The  degrees-of -f reedom  include  vertical  and  horizontal 
translation  (no  forward  translation),  roll,  pitch,  and  yaw.  A simpli- 
fied version  of  the  model  referred  to  as  TRIP  has  no  separate  sprung 
mass.  A more  complex  version  of  the  model,  referred  to  as  SPTRIP 
includes  separate  sprung  mass  roll,  lateral  translation  and  heave 
degrees-of -freedom.  Additional  model  components  simulate  tire  sliding 
and  normal  forces,  suspension  springs,  a lateral  sprung  mass  bushing 
force,  and  a mode  that  simulates  the  reaction  force  from  the  curb  due  to 
vehicle  impact.  For  a detailed  explanation  and  listing  of  the  equations 
and  modes  the  reader  is  referred  to  Appendix  A. 

B.  INPUT  AND  OUTPUT  FEATURES 

The  program  is  designed  for  user  ease.  The  vehicle  characteristics 
(track  width,  sprung  and  unsprung  masses,  tire  damping,  etc.)  are 
located  in  a file,  which  is  then  read  into  the  program  (refer  to  Appen- 
dix B).  The  user  must  only  enter  the  data  file  name  and  the  vehicles 
initial  dynamic  conditions.  The  program  allows  the  user  to  choose 
between  either  iterating  to  the  tip  over  velocity  (within  ±0.1  ft/sec) 
or  running  one  set  of  dynamic  conditions  and  seeing  the  outcome.  The 
program  also  allows  the  user  to  input  the  dynamic  variables  that  will  be 
plotted.  Any  combination  of  variables  for  any  set  of  conditions  can  be 
plotted.  The  user  also  must  choose  the  time  span  that  will  be  run  and 
the  iteration  time  increment.  Finally  the  user  is  allowed  to  input  one 
line  of  comments  which  will  be  included  in  the  printed  output.  The  pre- 
ceeding  input  features  will  be  further  discussed  in  Section  III-B. 

The  output  consists  of  transient  response  plots  of  the  four  varia- 
bles chosen  during  the  input  process.  These  plots  are  labeled  complete 


TR-1216-1 


6 


with  units  and  scales  which  will  appear  on  the  screen.  The  program 
automatically  prints  out  (to  a printer)  the  initial  and  impact  condi- 
tions as  well  as  the  line  of  comments. 

C.  TEST  CASES  FOR  SPTRIP 

(MODEL  INCLUDING  SPRUNG  MASS) 

This  program  is  for  the  simulation  of  a dynamic  vehicle  model, 
therefore  elementary  test  cases  were  run  and  evaluated  to  determine  if 
the  dynamic  equations  are  behaving  properly.  Twenty-six  different  test 
cases  were  run  to  establish  that  all  degrees-of -freedom  were  in  working 
order.  ine  ensuing  test  cases  were  run  using  vehicle  parameters  for  a 
1975  Chevy  Nova.  The  values  for  these  parameters  can  be  found  in  the 
back  of  Appendix  B.  Definitions  for  the  plotting  variables  are  given  in 
Section  III. 

1.  Vertical  Translation  (Heave) 

The  first  test  case  involves  a car  being  dropped  from  above  the 
ground  (none  of  the  tires  are  in  contact  with  the  ground).  The  purpose 
of  this  est  was  to  show  tnat  both  the  tires  and  suspension  springs  damp 
out  to  constant  values  over  a period  of  time.  Figure  4 shows  that  over 
a short  period  of  time  both  the  center  of  gravity  heights  and  velocities 
for  each  mass  do  indeed  damp  out.  Note  that  the  sprung  mass  has  a natu- 
ral frequency  on  the  order  of  2 Hz  (i.e.,  ZSD). 

2.  Lateral  Translation 

Cases  2 and  3 (Fig.  5)  show  the  car  sliding  both  to  the  left  and  the 
right.  The  purpose  for  this  was  to  show  that  the  car  slows  to  a stop 
due  to  sliding  friction  and  comes  to  rest  in  either  direction.  These 
plots  show  that  due  to  car  sliding,  small  roll  angles  appear,  as  was 
expected  due  to  the  lateral  deceleration  induced  by  the  lateral  tire/ 
road  friction  force. 
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Figure  4.  Test  Case  for  Sprung  and  Unsprung 
Mass  Vertical  Degrees-of-Freedom 
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COMMENTS  : TIRE  AND  SUSPENSION  DAMPING  (SPRUNG)  3/26/03 

THE  INITIAL  CONDITIONS  WERE: 

ZU  =-1.5  FT  ZUD  “ O FT/SEC  PHIU  = 0 DEG  PHIUD  = O'  DEG/SEC 

PHIS  ° O DEG  PHISD  = O DEG /SEC  THE  = O DEG  THED  = O DEG /SEC 

YO  - 10  FT  YUD  = O FT /SEC  YAWR  = 0 DEG /SEC  PS I - O DEG 
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Figure  5.  Test  Cases  for  the  Lateral  Slide 
Degree-of -Freedom 
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3.  Roll 


Cases  4-7  (Fig.  6)  simulate  the  car's  response  in  the  roll  degree- 
of-freedom.  In  all  the  cases  tested  a given  roll  angle  produced  a con- 
vergent roll  reaction  in  both  the  sprung  and  unsprung  masses.  In  the 
case  of  the  sprung  mass,  it  was  tilted  to  a 5 deg  angle  and  released 
while  the  unsprung  mass  angle  remained  zero.  Figures  6a  and  6b  show 
that  the  resulting  motion  of  the  system  converged  to  zero  for  both  posi- 
tive and  negative  initial  angles,  with  a natural  frequency  on  the  order 
of  2 Hz.  For  the  case  involving  the  unsprung  mass  the  car  was  tilted  to 
an  angle,  with  either  the  right  or  left  wheels  remaining  in  contact  with 
the  ground,  and  released.  Figures  6c  and  6d  show  that  like  the  sprung 
mass,  the  system  damps  out  for  both  positive  and  negative  initial 
angles. 

4.  Pitch 

Cases  8 and  9 (Fig.  7)  represent  the  pitching  degree-of -freedom. 
For  these  cases  either  the  front  or  the  rear  tires  are  lifted  off  the 
ground  (while  the  others  remain  in  contact)  and  are  released.  The 
resulting  motion  is  convergent  and  is  shown  in  Fig.  7a  for  positive  ini- 
tial pitch  angle,  and  in  Fig.  7b  for  a negative  angle. 

5.  Yaw  (Heading  Angle) 

Cases  10-13  are  a compilation  of  tests  involving  the  yawing  motion 
of  the  car.  In  Cases  10  and  11  (Figs.  8a  and  8b)  the  car  was  given  only 
an  initial  yaw  rate.  The  purpose  for  these  cases  was  to  show  that  the 
yaw  rate  decreased  to  zero  while  the  yaw  angle  approaches  a constant 
value.  This  test  was  run  for  both  positive  and  negative  yaw  rates,  and 
in  both  cases  the  results  were  found  to  be  dynamically  consistent  with 
expected  real  world  behavior.  The  other  yawing  condition  considered  was 
for  a car  sliding  at  an  initial  yaw  angle.  When  the  car  came  to  a stop 
the  yaw  angle  should  have  remained  the  same.  As  can  be  seen  in  Figs.  8c 
and  8d  this  was  indeed  the  case  for  both  positive  and  negative  yaw 
angles. 
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a)  Initial  Positive  Sprung  Mass  Roll  Angie 
Perturbation 
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b) Initial  Negative  Sprung  Mass  Roll  Angie 
Perturbation 
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Figure  6.  Test  Cases  for  Both  Unsprung 
and  Sprung  Mass  Roll  Perturbations 
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c)  Initial  Positive  Unsprung  Mass  Roll  Angle 
Perturbation 
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d ) Initial  Negative  Unsprung  Mass  Roll  Angle 
Perturbation 
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Figure  6.  (Concluded) 
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COMMENTS  l INITIAL  ROLL  ANGLE  (STRUNG  MASS)  3/26/B3  COMMENTS  : INITIAL  ROLL  ANGLE  (SPRUNG  MASS)  3/26/05 

THE  INITIAL  CONDITIONS  WERE  I THE  INITIAL  CONDITIONS  WEREt 

ZU  —1.3  FT  ZUD  * 0 FT/SEC  PHIU  = 10  DEG  PHIUD  = 0 DEG/SEC  ZU  — »-5  ET  ZUD  = O FT/SEC  PHIU  *-10  DEG  PHIUD  « O DEG/SEC 

PHIS  * 10  DEG  PHISD  “ 0 DEG/SEC  THE  =>  0 DEG  THED  » O DEG/SEC  PHIS  *-10  DEG  PHISD  * 0 DEG/SEC  THE  = 0 DE6  THED  = 0 DEG/SEC 

TO  * 10  FT  YUD  » 0 FT/SEC  YAWP  - 0 DEG/SEC  PS  I = 0 DEG  YO  = 10  FT  YUD  = O FT/SEC  YAWR  = 0 DEG/SEC  PS l “0  DEG 
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a)  Initial  Positive  Pitch  Angie 
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b)  Initial  Negative  Pitch  Angle 
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Figure  7.  Test  Cases  for  Pitch 
Degree- of- Freedom 
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Figure  8.  Test  Cases  for  the  Yaw 
Degree-of -Freedom 
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d)  Initial  Negative  Yaw  Angle 
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Figure  8.  (Concluded) 


TR-1216-1 


15 


6.  Curb  Impact 


Cases  14-17  (Fig.  9)  represent  initial  conditions  set  up  so  that  the 
car,  traveling  at  some  velocity  and  some  initial  heading  angle,  impacts 
the  curb.  Cases  14  and  15  were  run  with  a non-symmetric  car  (distance 
from  the  center  of  gravity  to  the  front  and  rear  axles  was  not  the 
same),  whereas  Cases  16  and  17  are  for  a symmetric  car  (distance  from 
the  center  of  gravity  to  the  front  and  rear  axles  was  equal).  For  the 
non-symme trical  case  when  the  car  impacts  the  curb  from  both  positive 
and  negative  heading  angles  (same  magnitude  in  each  case)  there  will  be 
some  minor  differences  due  to  the  non-symme trical  conditions.  This  can 
be  seen  in  Figs.  9a  and  9b.  For  symmetrical  cases  the  vehicle  response 
should  be  the  same  regardless  of  which  direction  the  heading  angle  acts. 
This  can  be  seen  in  Figs.  9c  and  9d. 

In  Cases  18  and  19  (Fig.  10)  the  car  slides  into  the  curb  with  only 
an  initial  velocity.  Figure  10a  represents  a case  were  the  car  did  not 
roll  over.  This  indicates  that  the  impact  velocity  was  not  sufficient 
enough  to  tip  the  car,  and  a higher  velocity  is  needed.  Figure  10b  once 
again  represents  a case  where  the  car  did  not  tip  over.  However, 
Fig.  10b  was  obtained  using  the  program's  iteration  method  to  determine 
the  tip  over  velocity.  Near  the  top  of  Fig.  10b  is  the  impact  velocity, 
this  velocity  (which  is  slightly  greater  than  that  shown  in  Fig.  10a)  is 
the  maximum  velocity  at  which  the  car  may  impact  the  curb  and  not  roll- 
over. This  brings  us  to  Fig.  10c  which  represents  the  car  tipping 
over.  The  impact  velocity  for  this  case  is  slightly  greater  than  that 
from  Fig.  10b,  yet  the  car  did  roll  over.  This  shows  that  the  iteration 
method  works  to  approximately  ±0.1  ft/sec.  The  tip  over  velocities 
shown  in  Fig.  10  are  in  the  range  of  that  expected. 

D.  TEST  CASES  FOR  TRIP  (MODEL  WITHOUT 

SEPARATE  SPRUNG  MASS) 

Before  the  final  version  of  the  program  was  developed  and  tested  a 
less  complex  car  model  was  developed.  This  model  did  not  include  the 
sprung  mass  degree-of -f reedom.  Instead  the  earlier  program  was  designed 
using  only  the  unsprung  mass  with  a center  of  gravity  height  equivalent 


TR-1216-1 


16 


a)  Non-Symmetric  (yz  plane)  With  Positive 
Heading  Angie 
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b)  Non-Symmetric  (yz  plane )W it h Negative 
Heading  Angie 
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Figure  9.  Test  Case  for  a Car  Impacting 
Curb  with  a Heading  Angle 
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COMMENTS  i SLIDING  INTO  CURB  WITH  HEADING  (SPRUNG  MASS)  3/26/B5  COMMENTS  : SLIDING  INTO  CURB  WITH  HEADING  (SPRUNG  MASS)  3/26/65 

THE  INITIAL  CONDITIONS  WERE  * THE  INITIAL  CONDITIONS  WERE: 

ZU  =-1.11  FT  ZUD  =-  0 FT/SEC  PHIU  = 0 DEG  PHIUD  = O DEG/SEC  zu  =-1-11  FT  ZUD  = O FT/SEC  PHIU  » 0 DEG  PHIUD  = 0 DEG/SEC 

PHIS  = 0 DEG  PH1SD  *=  0 DEG/SEC  THE  = O DEG  THED  = 0 DEG/SEC  PHIS  = 0 DEG  PHI5D  = 0 DEG/SEC  THE  = O DEG  THED  = 0 DEG/SEC 

VO  = 10  FT  YUD  = 22  FT/SEC  YAWR  = 0 DEG/SEC  PSI  = 10  DEG  TO  - 10  FT  YUD  = 22  FT/SEC  YAWR  = 0 DEG/SEC  PSI  =-10  DEG 


c ) Symmetric  ( yz  plane)  With  Positive 
Heading  Angie 
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d)  Symmetric  (yz  plane  ) With  Negative 
Heading  Angle 
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Figure  9.  (Concluded) 
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THE  INITIAL  CONDITIONS  WERE:  THE  INITIAL  CONDITIONS  WERE: 

ZU  =-1.11  FT  ZUD  = 0 FT/SEC  F'HIU  = 0 DES  FHIUD  = 0 DEG/SEC  zu  t_i.u  py  ZUD  = 0 FT/SEC  PHIIJ  = 0 DEG  PHIUD  = 0 DEG/SEC 

PHIS  - 0 DEG  PHISD  = 0 DEG/SCC  THE  = 0 DEG  THED  = 0 DEG/SEC  PHIS  = 0 DEG  F'HISD  = 0 DEG/SEC  THE  = 0 DEG  THED  = 0 DEG/SEC 

YT.  «=  In  FT  YUD  = 22  FT/SEC  YAWR  = 0 DEG/SEC  PSI  = 10  DEG  YO  = 10  FT  YUD  = 22  FT/SEC  YAWR  = 0 DEG/SEC  PS  I =-10  DEG 
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b)  Maximum  Case  Where  Car  Does  Not 
Tip  Over 
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Figure 


10.  Test  Case  for  Car  with  No 
Heading  Impacting  Curb 
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c)Car  Tip  Over 
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Figure  10.  (Concluded) 


to  that  of  the  sprung  and  unsprung  masses  together.  Once  the  sprung 
mass  program  was  finally  developed  it  was  useful  to  compare  the  two  pro- 
grams. Figures  11-13  represents  the  final  six  test  cases  that  were  run. 
Figures  11a  and  lib  show  that  the  tip  over  angle  for  the  unsprung  mass 
program  (TRIP)  was  between  45  and  46  deg,  which  was  expected.  For  the 
sprung  mass  program  (SPTRIP)  the  suspension  springs  were  stiffened  and 
Figs.  12a  and  12b  show  a tip  over  angle  between  49  and  50  deg.  The  dif- 
ference in  tip  over  angle  is  attributed  to  the  fact  that  although  the 
suspension  springs  were  stiffened,  it  was  impossible  to  make  them  com- 
pletely rigid  and  therefore  they  deflected  and  caused  varying  forces. 
It  should  be  noted  here  that  the  tip  over  angle  from  Figs.  12a  and  12b 
is  not  the  same  as  that  in  Fig.  10  due  to  the  fact  that  when  the  car 
impacts  with  the  curb  the  axle  is  crushed  and  the  half  track  width  is 
shortened,  therefore  making  tip  over  easier.  In  the  cases  shown  in 
Figs.  11-12  the  left  wheels  were  lifted  off  the  ground  and  released, 
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a)  No  Tip  Over 
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Figure  11.  Test  Case  for  the  Static  Tip  Over 
Condition  Using  the  Unsprung 
Mass  (TRIP)  Program 
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COMMENTS  : COMMENTS  : 

THE  INITIAL  CONDITIONS  WERE:  THE  > CONDITIONS  WERE: 

2L.I  =-3,44  FI  ZUD  = O FT/SEC  FHIU  = 45  DEG  ZU  —i,  .4  FT  ZUD  = O FT/SEC  PHIU  = 46  DEG 

PHIUD  “ 0 DEG/SEC  VO  » 10  FT  VUD  = 0 FT/SEC  PHILJD  « O DEG/SEC  YO  - 10  FT  YUD  = O FT/SEC 

THE  = O DEG  THED  = O DEG/SEC  YAWR  « O DEG/SEC  FSI  = O DEG  THE  = O DEG  THED  = 0 DEG/SEC  YAWR  = 0 DED/SEC  F'SI  = 0 DEG 
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b)  Vehicle  Tip  Over 
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Figure  12.  Test  Case  for  the  Static  Tip  Over 
Condition  Using  the  Sprung  Mass 
(SPTRIP)  Program 
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COMMENTS  : COMMENTS  : 

THE  INITIAL  CONDITIONS  WERE:  THE  INITIAL  CONDITIONS  WERE: 

ZU  =-2.65  FT  ZUD  = O FT/SEC  PHIIJ  = 49  DEG  PHIIJD  = O DEG/SEC  ZU  =-2.65  FT  ZUD  = O FT/SEC  PHIU  = 50  DEG  PH1UD  = O DEG/SEC 

PHIS  = 49  DEG  PHISD  = O DEG/SEC  THE  « 0 DEG  THED  = O DEG/SEC  PHIS  = 50  DEG  PHISD  = O DEG/SEC  THE  = O DEG  THED  = O DEG/SEC 

YO  = 10  FT  YUD  » O FT/SEC  YAWR  » 0 DEG/SEC  FSI  = O DEG  YO  = 10  FT  YUD  = O FT/SEC  YAWR  = 0 DEG/SEC  PSI  » 0 DEG 


a)  Car  Sliding  into  Curb  Using  TRIP 


b)  Car  Sliding  into  Curb  Using  SPTRiP 
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Figure  13.  Comparison  Between  the  Unsprung 
Mass  (TRIP)  and  Sprung  Mass 
(SPTRIP)  Programs 
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therefore  the  axles  did  not  crush  and  tip  over  was  more  difficult.  The 
final  two  cases  (Figs.  13a  and  13b)  are  provided  to  show  the  similarity 
between  the  two  programs  when  the  car  impacts  the  curb.  Both  Figs.  13a 
and  13b  reach  the  same  roll  angle  but  with  different  impact  velocities. 
It  can  be  seen  that  for  the  sprung  mass  program  (Fig.  13b)  a greater 
velocity  was  needed  to  produce  the  same  roll  angle  as  the  unsprung  mass 
program.  These  results  show  that  the  two  programs  are  similar,  and 
adding  the  sprung  mass  degrees-of-f reedom  did  not  adversly  affect  the 
model. 
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SECTION  III 


USER  INSTRUCTIONS 


A.  INSTALLATION  AND  OPERATION 


This  program  was  designed  to  operate  on  an  IBM®  personal  computer 
(or  any  compatible  system).  The  disk  is  set  up  for  easy  operation  so 
that  even  with  little  or  no  computer  experience  the  program  is  easily 
accessible. 


By  placing  the  program  diskette  in  drive  A (the  left  disk  drive)  and 
either  turning  the  computer  on  or  booting  the  system,  i.e.,  (simulta- 
neously depressing  the  control  Ctrl  , alternate  | Alt  | , and  delete 
keys)  the  computer  will  startup,  load  the  system,  and  install  the 


Del 


graphics  routine.  From  here  the  user  must  decide  which  form  of  the  pro- 
gram will  be  used. 


The  original  disk  contains  a compiled  version  of  both  the  unsprung 
mass  program  (TRIP)  and  the  sprung  mass  program  (SPTRIP).  These  com- 
piled versions  will  only  work  with  computers  that  have  an  8087  math 
coprocessor  chip.  For  this  reason  the  Basic  source  code  has  also  been 
included  on  the  disk  so  that  users  may  compile  the  program  to  suit  their 
needs.  The  two  compiled  versions  are  filed  under  the  names: 


TRIP8  7.EXE 
SPTRIP87.EXE 


When  the  prompt  appears  (A>)  typing  in  the  appropriate  file  name  and 
entering  it  will  execute  the  program. 


Finally  the  capital  lock 


Caps 

Lock 


key  must  be  activated  so  that  capi- 
tal letters  appear  when  keys  are  depressed.  If  not,  when  the  program  is 


running  it  will  not  execute  properly. 
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B.  USER  INPUTS  AND  PLOTTING  OUTPUT 


The  initial  conditions  that  the  user  must  enter  into  the  computer 
during  program  operation  will  now  be  discussed.  The  following  pages 
show  the  actual  input  questions  (underlined  terms  indicate  user  inputs), 
that  will  be  asked  during  program  operation.  These  will  be  used  as  a 
reference  when  discussing  each  individually.  The  following  discussion 
will  focus  on  the  sprung  mass  program  (SPTRIR).  The  inputs  for  the 
unsprung  mass  program  (TRIP)  are  similar  with  the  exception  that  they  do 
not  contain  the  sprung  mass  degree  of  freedom  inputs  and  options. 

User/computer  dialog  is  shown  explicitly  in  Table  1. 

The  first  prompt  in  Table  1 instructs  the  user  to  enter  the  name  of 
the  file  that  is  to  be  used  during  program  execution.  The  user  is 
reminded  to  use  both  a disk  drive  specifier  (A:,  B:,  etc.)  and  extension 
(.75,  .BAS,  etc.)  if  they  apply. 

The  second  Table  1 prompt  allows  the  user  to  either  iterate  to  the 
vehicles  tip  over  condition  using  velocity  as  the  variable,  or  to  look 
at  one  particular  set  of  conditions  at  a time. 


The  third  Table  1 prompt  is  dependent  on  the  form  of  the  previous 
input.  As  was  mentioned  earlier  if  the  inputs  are  not  in  capital  let- 
ters the  program  will  not  function  properly.  If  the  answer  to  question 
2 was  input  using  capital  letters  the  program  will  bypass  this  section 
and  proceed.  If  the  answer  was  input  using  lower  case  letters  this  sec- 


tion  will  prompt  the  user  to  first  activate  the  capital  lock 
and  then  reenter  the  answer. 


Caps 

Lock 


key 


The  fourth  Table  1 prompt  designates  which  four  variables  the  user 
would  like  to  see  on  the  plot,  and  what  the  scaling  for  each  will  be.  A 
list  of  variables  will  appear  with  a short  explanation  of  what  each 
variable  represents.  From  the  list,  four  variables  are  choosen  and 
entered  one  at  a time  with  their  desired  scale  (e.g.  , PHIS, 40  or  YUD, 20 , 
etc. ) . 
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TABLE  1.  USER/COMPUTER  DIALOG  (UNDERLINED  TERMS 
INDICATE  USER  INPUTS) 


1.  ENTER  THE  DESIRED  FILE  NAME 

(INCLUDE  DRIVE  SPECIFIER  IF  NEEDED) 
? A: SPNOVA. 7 5 ? 


2.  WOULD  YOU  LIKE  THE  PROGRAM  TO  ITERATE  TO 
THE  TIP  OVER  VELOCITY  (Y  OR  N)? 

? n ? 


3.  PLEASE,  DEPRESS  THE  CAPS  LOCK  KEY  (RIGHT  OF  THE  SPACE  BAR) 
ONE  TIME  AND  THEN  REENTER  YOUR  REPLY  TO  THE  ABOVE 
? N 


4.  PLOTTING  VARIABLE  SELECTION: 


PHIS 

- 

THE 

SPRUNG  MASS  ROLL  ANGLE 

PHISD 

- 

THE 

SPRUNG  MASS  ROLL  RATE 

PHISDD 

- 

THE 

SPRUNG  MASS  ROLL  ACCELERATION 

PHIU 

- 

THE 

UNSPRUNG  MASS  ROLL  ANGLE 

PHIUD 

- 

THE 

UNSPRUNG  MASS  ROLL  RATE 

PHIUDD 

- 

THE 

UNr  'RUNG  MASS  ROLL  ACCELERATION 

PSI 

- 

THE 

YAk  ANGLE  AS  THE  CAR  APPROACHES  THE  CURB 

THE 

- 

THE 

CAR'S  PITCH  ANGLE 

THED 

- 

THE 

CAR'S  PITCH  RATE 

THEDD 

- 

THE 

CAR'S  PITCH  ACCELERATION 

YAWR 

- 

THE 

RATE  OF  YAW  AS  THE  CAR  APPROACHES  THE  CURB 

YAWRD 

- 

THE 

YAW  ACCELERATION  AS  THE  CAR  APPROACHES  THE  CURB 

YO-YU 

- 

THE 

c.g.'s  DISTANCE  FROM  CURB 

YU 

- 

THE 

c.g.'s  POSITION  WITH  RESPECT  TO  IT'S  INITIAL  CON 

YUD 

- 

THE 

c.g.'s  LATERAL  VELOCITY 

YUDD 

- 

THE 

c.g.'s  LATERAL  ACCELERATION 

ZU 

- 

THE 

c.g.'s  VERTICAL  POSITION  WITH  RESPECT  TO  THE  GROUND 

ZUD 

- 

THE 

c.g.'s  VERTICAL  VELOCITY 

ZUDD 

- 

THE 

c.g.'s  VERTICAL  ACCELERATION 

DELZU 

— 

THE 

CHANGE  IN  ZU  FROM  ITS  STATIC  POSITION 

ENTER  4 

VARIABLES  AND  THEIR  SCALES  ONE  AT  A TIME  (i.e.,  PHIS, 40) 

? PHIS, 40  7 
? PHIU, 40 7 
? YUD , 20  ? 

? YO-YU, 10  y 
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TABLE  1.  (CONCLUDED) 


5.  INITIAL  CONDITION  SPECIFICATION  FOR  VARIABLES: 

ENTER  : ZU  (FT),  ZUD  (FT/SEC),  PHIU  (DEG),  PHIUD  (DEG/SEC) 

? -1.11,0,0,0  ? 


ENTER  : PHIS  (DEG),  PHISD  (DEG/SEC),  THE  (DEG),  THED  (DEG/SEC) 

? 0,0, 0,0  p 


ENTER  : YO  (FT),  YUD  (FT/SEC),  PSI  (DEG),  YAWR  (DEG/SEC) 
? 10 , 22  ; 0 , 0 ;? 


6.  COMPUTATION  TIME  INCREMENT  AND  TIME  SPAN: 

ENTER  THE  DESIRED  TIME  INCREMENT  (<=.0025) 
? .0025  y 


ENTER  THE  DESIRED  TIME  LIMIT  (INTEGERS  ONLY) 
? 2 * 


ENTER  ONE  LINE  OF  COMMENTS  (PLEASE,  NO  COMMAS)? 

? ANY  LINE  OF  COMMENTS  CAN  GO  IN  HERE  BUT  NO  COMMAS  y 
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The  fifth  set  of  prompts  requests  specification  of  the  vehicle's 
initial  conditions.  There  are  several  important  things  to  remember  when 
entering  the  initial  conditions.  The  first  is  that  Zu  (the  unsprung 
mass  center  of  gravity  vertical  distance,  above  the  ground)  should  be 
negative  and  close  to  the  trim  value  for  the  vehicle,  and  that  it 
depends  on  the  initial  pitch  and  roll  angles.  Secondly  YQ  (the  initial 
perpendicular  distance  from  the  center  of  gravity  to  the  curb)  must 
start  at  least  the  distance  of  one-half  the  track  width  from  the  curb. 
Otherwise  the  vehicle  will  be  past  the  curb  at  the  start. 

Due  to  the  complexity  of  the  program  and  the  different  dynamic  modes 
it  contains  it  was  found  to  run  best  when  time  increments  of  0.0025  and 
under  are  used.  The  rest  of  the  initial  vehicle  conditions  are  user 
dependent  and  the  only  restriction  is  that  they  are  entered  in  the  spec- 
ified units.  The  desired  time  limit  refers  to  the  maximum  amount  of 
simulated  and  plotted  time  that  the  user  would  like.  This  input  must  be 
input  in  seconds  and  only  integer  values  should  be  used. 


The  final  input  question  allows  the  user  to  input  a line  of  comments 
that  will  appear  on  the  plots.  Due  to  the  structure  of  the  program, 
however,  no  commas  may  be  used. 


Once  all  of  the  inputs  have  been  entered  a grid  will  appear  on  the 

screen  and  the  curves  will  be  plotted.  If  the  program  is  iterating, 

each  iteration  will  produce  a new  set  of  time  histories  for  the  user  to 

see.  When  the  program  is  through  running  it  will  beep  to  alert  the  user 

that  it  is  done.  The  program  will  automatically  print  out  the  line  of 

comments,  the  initial  conditions,  and  the  conditions  upon  impact.  At 

this  point  the  user  can  either  get  a copy  of  the  variable  time  histories 

or  change  the  dynamic  inputs  or  plotting  variables.  To  obtain  a copy  of 

the  time  histories  the  user  must  simultaneously  depress  the  shift  ^ 

Pr t Sc 

and  print  screen  ^ keys.  After  the  print  is  done  or  before 
printing,  depressing  any  key  will  return  the  user  to  the  next  section  of 
the  program. 
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C.  REINITIALIZATION 


This  section  of  the  program  allows  the  user  to  change  any  of  the 
initial  dynamic  inputs  (Table  1,  prompt  5),  and  (or)  the  plotting  vari- 
ables (Table  1,  prompt  4). 

The  program  will  first  prompt  the  user  and  ask  if  any  of  the  initial 
dynamic  variables  are  to  be  changed,  Table  2,  prompt  1.  If  the  user’s 
response  is  Y_  (Yes),  the  program  will  list  the  variables  that  may  be 
changed  and  then  prompt  (Table  2,  prompt  2)  the  user  to  enter  the  vari- 
able name,  followed  by  a comma,  and  then  the  new  value  for  that  variable 
(e.g.,  YUD,  30).  After  the  change  is  entered  the  program  will  prompt 
the  user  for  another  change.  This  process  continues  until  all  changes 
are  made.  At  this  point,  (or  any  other  time)  when  the  program  prompts 
for  the  next  change,  the  user  may  continue  by  entering  a Q_  , followed  by 
a comma  (e.g.,  Q_,).  The  program  will  now  proceed  to  the  section  that 
allows  the  user  to  change  the  plotting  variables.  It  should  be  noted 
that  entering  the  N_  (No)  option  when  the  program  asks  if  any  dynamic 
variables  are  to  be  changed  (Table  2,  prompt  2),  will  automatically  take 
the  user  to  the  section  that  allows  the  plotting  variables  to  be 
changed. 

The  program  will  now  prompt  the  user  and  ask  if  any  of  the  plotting 
variables  will  be  changed,  (Table  2,  prompt  3).  If  the  user's  response 
is  Y_  (Yes),  the  program  will  list  all  of  the  possible  plotting  vari- 
ables. The  user  may  now  change  the  plotting  variables,  (Table  2,  prompt 
4)  in  the  same  manner  as  was  shown  in  Table  1,  prompt  4.  It  should  be 
noted  that  if  the  user  decides  to  change  any  plotting  variable,  four 
variables  must  be  input  not  just  one.  After  these  four  variables  and 
their  scales  are  input  the  program  will  automatically  plot  the  new  time 
his  tories. 

Finally,  a N_  (No)  response  to  prompt  3 will  cause  one  of  two  events 
to  occur.  If  prompt  1 is  yes  and  prompt  3 is  no,  the  program  will  plot 
the  new  time  histories  using  the  previous  plotting  variables.  If  prompt 
1 is  no  and  prompt  3 is  no  the  program  will  cease  execution. 
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TABLE  2.  CHANGING  INITIAL  INPUTS  (UNDERLINED  TERMS 
INDICATE  USER  INPUTS) 


1.  WOULD  YOU  LIKE  TO  CHANGE  ANY  OF  THE  INITIAL  DYNAMIC 
VARIABLES  (Y  OR  N)? 

?Y  ;? 


2.  YOU  MAY  CHANGE  ANY  OF  THE  VARIABLES  LISTED  BELOW: 


zu 

ZUD 

PHIU 

PHIUD 

PHIS 

PHISD 

THE 

THED 

YO 

YUD 

PSI 

YAWR 

ENTER  THE  VARIABLE  THAT  YOU  WOULD  LIKE  CHANGED,  FOLLOWED  BY  THE  NEW 
NUMERIC  VALUE  FOR  THE  VARIABLE  (i.e.,  YUD,  25),  OR  ENTER  (Q,)  TO 
CONTINUE 
?YUD , 30  ;? 


ENTER  ANOTHER  VARIABLE  TO  CHANGE  OR  ( Q , ) TO  CONTINUE 

?Q,  ^ 


3.  WOULD  YOU  LIKE  TO  CHANGE  THE  PLOTTING  VARIABLES  (Y  OR  N)? 
?Y  ^ 


4. 


PHIS 

- 

THE 

PHISD 

- 

THE 

PHISDD 

- 

THE 

PHIU 

- 

THE 

PHIUD 

- 

THE 

PHIUDD 

- 

THE 

PSI 

- 

THE 

THE 

- 

THE 

THED 

- 

THE 

SPRUNG  MASS  ROLL  ANGLE 

SPRUNG  MASS  ROLL  RATE 

SPRUNG  MASS  ROLL  ACCELERATION 

UNSPRUNG  MASS  ROLL  ANGLE 

UNSPRUNG  MASS  ROLL  RATE 

UNSPRUNG  MASS  ROLL  ACCELERATION 

YAW  ANGLE  AS  THE  CAR  APPROACHES  THE  CURB 

CAR'S  PITCH  ANGLE 

CAR'S  PITCH  RATE 
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TABLE  2.  (CONCLUDED) 


THEDD 

- 

THE 

YAWR 

- 

THE 

YAWRD 

- 

THE 

YO-YU 

- 

THE 

YU 

- 

THE 

YUD 

- 

THE 

YUDD 

- 

THE 

ZU 

- 

THE 

ZUD 

- 

THE 

ZUDD 

- 

THE 

DELZU 

- 

THE 

CAR'S  PITCH  ACCELERATION 

RATE  OF  YAW  AS  THE  CAR  APPROACHES  THE  CURB 
YAW  ACCELERATION  AS  THE  CAR  APPROACHES  THE  CURB 
c.g.'s  DISTANCE  FROM  CURB 

c.g.'s  POSITION  WITH  RESPECT  TO  IT'S  INITIAL  CON. 
c.g.'s  LATERAL  VELOCITY 
c.g.'s  LATERAL  ACCELERATION 

c.g.'s  VERTICAL  POSITION  WITH  RESPECT  TO  THE  GROUND 

c.g.'s  VERTICAL  VELOCITY 

c.g.'s  VERTICAL  ACCELERATION 

CHANGE  IN  ZU  FROM  ITS  STATIC  POSITION 


ENTER  4 VARIABLES  AND  THEIR  SCALES  ONE  AT  A TIME  (i.e.,  PHIS, 40) 
? PHIU , 40  y 
? PHIS, 40  7 
? YUD , 20  ? 

? YO-YU, 10  y 
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APPENDIX  A 


EQUATIONS  OF  MDTION  FOR  TRIPPED  ROLLOVER  ANALYSIS 

1.0  GENERAL  APPROACH  AND  OVERVIEW 

The  vehicle  degrees  of  freedom  include  for  the  sprung  and  unsprung 
mass:  roll  <t>u),  yaw  (^)»  and  pitch  (0)  angles,  lateral  position 

(y),  and  vertical  position  (z).  Constraints  on  the  sprung  relative  to 
the  unsprung  roll  angle  include  a torsional  spring  gradient  with  damping 
a"  bump  stops.  The  tires  include  vertical  compliance  and  lateral 
dynamic  friction.  The  curb  encounter  will  involve  a force/deflection 
characteristic.  The  force  and  moment  equations  will  be  developed  for 
both  the  sprung  and  unsprung  masses.  The  following  represents  the  vehi- 
cle model  used: 


TR-1216-1 


A-l 


The  axis  system  that  is  used  is  the  standard  aircraft  body  axis  conven- 
tion as  shown: 

+ x 


z 


+ 2 


+y 


•••* 


J U v 

> m $g$ 
f $ ft  H 


where  angles  are  measured  as  follows: 
For  pitch 


Front 

fY  Rear 

T7  /’f/’/'/’/'/'/’  2 /»/ 


For  roll 


For  yaw 
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1.1  SPRUNG  MASS  ROLL  AXIS 


By  applying  Newton's  second  law  (F  = ma)  to  the  following  free  body 
diagram,  both  force  equations  and  the  moment  equation  for  the  sprung 
mass  can  be  written. 


m m * t 

^ -.-V  -W.  JT-; 


:*S-  *V 


► y 

y 4>i  4>i4> 

z ->  defined  as  zero  when 
at  ground  level 


1.1.1  Lateral  Force  Equ  .ion 

“sYs  = (fSL  + fSR)  * sin  ~ F1  * cos  <j>s 

1.1.2  Vertical  Force  Equation 

- msg  + mszs  = - (Fsl  + Fsr)  • cos  (^  - Fi  • sin  4>s 

1.1.3  Moment  Equation 


Is^s  = (FSL  ~ FSR)  * s + Fi  • hp  1 
hFx  = (ys  " yu)  cos<|)s  + (Zu  - Zs)  cos<^  - (hr  - hu)  cos  ( ^ ) 

1.1.4  Suspension  Forces 

The  vehicle's  suspension  system  consists  of  both  left  and 
right  suspension  springs  and  bump  stops.  As  the  sprung  mass  and 
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unsprung  mass  roll  angles  change,  the  suspension  springs  expand  and 
compress  causing  suspension  force  changes.  The  bump  stops  prevent 


excessive 

tension  or 

compression  on  the  springs.  The  spring  force 

equations 

can  be  written  as: 

fsl  - K1 

• SD^  + B^  • SD^  + 

K.2  • BD-^  + B2  • BDl 

fSR  = K1 

• SD^  + Bj^  • SDr  + 

K2  • BDj^  + B2  • BD]^ 

Spring  Force  and 

Bump  Stop  Force  and 

Damp ing 

Damp  ing 

using: 

fs 

K2  / 

Forces  include  spring, 

Kl  J'' 

i 

1 * 

bump  stop  and  damping 

i 

1 

Spring 
i Deflection 

/ f 

i 

Bump 

Stop 

The  spring  deflections  can 

be  found  by  subtracting  the 

length  of 

the  spring 

at  any  instant  in 

time  from  the  unloaded  spring 

length. 

The  lengths 

of  the  springs  at  any  instant  in  time  can  be  found 

by  the  following: 
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Due  to  the  bump  stops,  Acjsg  can  not  reach  angles  much  greater  than  8°, 
therefore  it  can  be  reasonably  assumed  that  X is  perpendicular  to  the 
axle.  With  this  assumption  X can  be  written  as 

1/2 

X = [(zs  - zu)2  + (ys  - yu)2  ] + hu 


For  convenience  in  future  equations  the  following  notation 
will  be  employed. 


Az 

= zu 

" z s 

Az 

= zu 

- zs 

= Xs 

- yu 

Ay 

= ys 

- yu 

The  right  spring  length  can  be  written  as 


SR  = 7T~  ~ hs  “ a 

cos  A 4^ 


where 


a = s • tan  Acf>s 
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Therefore 


SR 


r 2 2 11/2 

[ tsz.1  + Ay‘  ] + hu 

cos  A<f>s  s 


s • tan  A4>s 


The  left  spring  length  can  now  be  written  as 
SL  = SR  + 2s  • tan  Acf>s 


Therefore 


SL 


[ Az^  + Ay^  ] + hu 

cos  Acfjg 


- hs  + s 


tan  A<j>s 


From  this  the  spring  deflection  terms  can  be  written  as 


SD^  = Spring  Length  - SL 


(A z » Az  4-  Ay 
cos  A<J>S  (Az^  + 


1/2 

A<(k  [s  + [(Az^  + Ay  2)  + hu]  sin  Acf>s  ] 


cos' 


and 

SDr  = Spring  Length  - SR 


SDr 


1/2 


(Az  » A?.  + Ay  * Ay ) [K  ^z2  + Ay2)  hul  s in  A<f>s  s ] 


cos  A^  (Az^  + Ay^) 


1/2 


cos  2 A<t£ 


The  bump  stop  deflections  are  found  directly  from  the  spring 
lengths  at  any  instant  in  time.  They  are  determined  by  subtracting  the 
spring  length  from  the  size  of  the  bump  stops.  If  the  spring  length  is 


A- 6 
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less  than  the  bump  stop  size,  then  the  axle  is  in  contact  with  the  bump 
stops.  The  size  of  the  bump  stops  is  determined  from  the  following: 


Bump  Stop  Size  = Static  Loaded  Spring  Length  - Static  Allowable  Bump 
Stop  Distance 


Therefore 


BDl  = Bump  Stop  Size  - SL 

bdl  = 0 

BDR  = Bump  Stop  Size  - SR 

BDr  = 0 


for  Bump  Stop  Size  > SL 
for  Bump  Stop  Size  < SL 
for  Bump  Stop  Size  > SR 
for  Bump  Stop  Size  < SR 


Since  the  bump  stop  sizes  are  essentially  constant  for  any  particular 
vehicle  it  follows  that 


BDl  = S*L 
BDr  = SR 


for  Bump  Stop  Size  > SL,  SR 


1.1.5  Lateral  Restraining  Force 

The  sprung  and  unsprung  masses  are  connected  together  by  a 
rubber  bushing,  therefore  causing  a lateral  force  between  the  two 
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masses. 


The  force  can  be  approximated  by  a spring  connection  at  this 


point. 


The  restraining  force  can  be  written 


= K • Def  + B • Def 


Using  the  above  illustration  Def  can  be  determined  as  follows: 


Def  = Zy’  - (hr  - hu)  sin 
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where 


Ay ' = 

[ Ay 

- Az 

• 

tan  <j^  ] 

cos 

4>s 

Ay ' = 

Ay 

• cos 

— A 7.  • 

sin 

4>s 

Therefore 

Def  = 

Ay 

• cos 

— \7.  • 

sin 

- (hr 

- hu)  sin  A^ 

and 

Def  = 

Ay 

• cos 

♦s 

• 

— Az  • 

sin 

‘t’s 

- k 

[Ay  • sin  <j)s  + 

- (hr  - hu)  A^  • cos  A<j^ 


1.2  UNSPRUNG  MASS  ROLL  AXIS 

By  applying  Newton's  second  law  to  the  following  free  body  diagram, 
both  force  equations  and  the  moment  equation  can  be  written  for  the 
unsprung  mass. 


A- 9 
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1.2.1  Lateral  Force  Equation 


Mu  Yu  = F1  cos  ~ (fSL  + fSR ) sin  4>S  “ fYR  “ fYL  " Fd 

1.2.2  Vertical  Force  Equation 

- mug  + muZu  = (FSl  + fSr)  cos  + F1  sin  - FZL  - FZR 


1.2.3  Moment  Equation 


i 

i 


i 


= <£s  <£u 
FZL  = FZLF  + FZLR 
FZR  = F ZRF  + FZRR 


When  one  or  both  tires  are  on  the  ground  the  moment  arm  for 
the  tire  sliding  forces  (FyR  ; FyR)  Fs  iust  simply  zu.  When  the  car 
comes  in  contact  with  the  curb  the  moment  arm  for  the  curb  forces  (F^) 
is  simply  zu  + Hcurb. 

Therefore 

IuJki  = ~ fSL  * s * cos  + Fsr  • s • cos  A<fe  + FZL  • y - FZR  • x 

“ fYL  * zu  “ fYR  * zu  + f1  * hrc  * cos  A<t> s 
" Fd  * <zu  + HCURB) 
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When  the  vehicle  begins  to  roll  the  moment  arms  for  the  vertical  tire 
forces  change  as  a function  of  the  tire  deflection  and  the  unsprung  mass 
roll  angle.  The  car's  track  width  (TRW)  is  divided  into  two  parts,  the 
left  (TRWL/2)  and  the  right  (TRWR/2).  This  is  due  to  the  fact  that  they 
may  be  shortened  due  to  structural  deformation  when  contacting  the  curb. 
This  all  can  be  seen  in  the  following: 


x = (TRWR/2  - a)  cos  4>u 
a = (hu  + ZTr)  tan  4^ 

Therefore 

x = (TRWR/2)  cos  4\i  - (hu  + ZTR)  sin  4>u 
y = TRW  • cos  4\j  - x 

Therefore 

y = (TRWL/2)  cos  4^  + (hu  + ZTR)  sin  4^ 
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Therefore 


Iu<ki  = (FSR  “ FSL)  * s * cos  A<fe  + fZL  [(TRWL/2)  cos  + (hu  + ZTR)  sin  (J^  ] 
” fZR  [(TRWR/2)  cos  <j^  - + ZTR)  sin  ^ ] - FYL  • zu  - • zu 

+ Fi  • \c  * cos  - Fd  * (zu  + HCURb) 

Caution:  The  TRW/2  term  may  be  shortened  due  to  structural  deformation 
due  to  contact  with  the  curb.  Therefore  in  this  equation  it  is  written 
for  the  particular  side  referenced  (TRWR/2  or  TRWL/2). 

1.3  YAW  AXIS  (SPRUNG  AND  UNSPRUNG) 


1.3.1  Moment  Equation 


!zr  - - (FYLF  + Fyrf  + FdF)  a 


+ (fYLR  + fYRR  + FdR)  b + ^s^Ix  “ 1y'> 
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1.4  PITCH  AXIS  (SPRUNG  AND  UNSPRUNG) 


dz 


Ix  ) r <^u 


gyroscopic 

moment 


1.4.1  Moment  Equation 


Iye  - (Iz  - Ix)  r • <tfe  - b(FZLR  + FZRr)  + a^FZLF  + FZRf) 


1.5  TIRE/WHEEL  FORCES 


1.5.1  Tire  Sliding  Forces 


For  the  sliding  friction  model.  The  side  force  is  a func- 
tion of:  normal  load,  coefficient  of  friction,  and  side  velocity,  and 
can  be  expressed  as 


fYLF  = b * 
FYRF  = y • 

fylr  = ^ * 

fyrr  = y * 

where  fy(v)  is  determined 


FZLF 

fZRF 

fZLR 

fZRR 

from: 


fy(yu  + a 

fy(yu  + a 

fy(yu  “ b 
fy(yu  “ b 


r - Roll  Term) 
r - Roll  Term) 
r - Roll  Term) 
r - Roll  Term) 


TR-1216-1 


A-l  3 


where, 

v = yu  +a*r-  Roll  Term 

'•  »•««#.  W-  - I 

or, 

v = Roll  Term 

..C  ' .>.«$  WV  .1. . / X&  ' 

■■  f*%  m "r#A'  f '-:s  ; figf  fe 

- m m ~ 'i  - > '%  m 

-•••  - & 


The  roll  term  is  the  velocity  caused  by  4^  as  the  vehicle  is  rolling. 
This  is  evaluated  as  follows: 


Roll  Term 


U<hu  + ZT)2  + (- 


TRW 


V? 

] cos 


$ 


where 


6 = 90  - a - 4^ 


and 


a = tan 


-1 
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1.5.2  Tire  Normal  Forces 


The  normal  forces  result  from  tire  deflections  determined  by 
the  distance  between  the  axle  and  the  road  surface  at  each  wheel.  The 
distance  depends  on  the  c.g.  height,  zu,  and  vehicle  pitch  and  roll 
attitude  (0,  <J>  ) . The  forces  include  both  spring  constant  and  damping 

terms . 

Tire  Deflections  (Compression  of  Tire) 

For  the  right  tires: 


Therefore 

ZTr  = — — - h - tan  <}>,,  + Pitch  Term 

R cos  u 2 ^ 

For  the  left  tires: 


ZTl  = ZTr  + TRW  • tan  ^ 
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Therefore 


-z„  TRW 

ZTt  = — - hn  + — — tan  ck,  + Pitch  Term 

^ cos  4^  u 2 u 

For  the  pitch  terms: 


For  the  left  and  right  front  tires  the  pitch  term  is 


a tan  0 


For  the  left  and  right  rear  tires  the  pitch  term  is 


- b tan  0 


Therefore  the  axle  heights  above  the  ground  are 


ZT 


LF 


ZT 


LF 


ZT 


RF 


Z*T 


RF 


"zu  i_  , TRW  ^ 

“ h..  + — * tan  <k,  + a • tan  0 

cos  (J)q  u 2 ^ 

~ zu  <Kj  (TRW/ 2 — zu  • sin  4>y)  + a • 0 

cos  (ki  cos^  (jjy  cos^  0 

— - h,,  - • tan  4.  + a • tan  0 

— zu  'Ki  (TRW/ 2 + zu  • sin  <J>y)  ^ a • 0 

cos  ^ cos^  4^  cos^  0 
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tan  4^  - b • tan  0 


ZT 


LR 


-z, 


- h„  + 


TRW 


cos  4^  u 2 


ZT 


LR 


4>u  (TRW/2  Zu  • sin  4>u)  b • 0 


cos  ‘tb  cos^  4^ 


cosz  0 


ZT 


-z, 


RR 


- h„  - 


TRW 


cos  4>u  u 2 


tan  4>u  ~ ^ * tan  9 


ZT 


-zu  4>u  (TRW/ 2 + zu  • sin  4^)  5.9 


RR 


cos  <ki  cos^  4^ 


cos/  0 


Using  the  axle  heights  and  the  plot  from  above,  the  force  for  each  tire 
can  be  written  as: 

FZT  - MZ<ZT>  + B2  • ZT 


where 


Kz  and  Bz  = 0 for  ZT  > Tire  Radius 


Fzt  = Kz  [(Tire  Radius)  - ZT^]  - Bz  • Z*T^ 
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1.5.3  Road  Surface  Lateral  Force  vs.  Displacement 
Characteristic  at  the  Discontinuity 


Road  surface  discontinuity  is  assumed  to  cause  both  elastic 
and  plastic  deformation  in  wheel/suspension  structure. 


yQ  = initial  conditions  for 
starting  distance  of 
c.g.  from  curb 

Body,  Frame,  Engine 


Contact 

With  Curb  TRWF  TRWR  Decreasing  > 
2 ’ 2 Values 
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yF  = y0  - yu  - a • sin  f - (TRWF/2)  cos 


y ? = - yu  - a • r • cos  ^ + (TRWF/2)  r • cos 


yR  = yQ  - yu  + b * sin  ^ “ (TRWR/2)  cos  ijj 


y ^ = - yu  + b • r • cos  ^ + (TRWR/2)  r • sin  \p 


FdF 


FdR 


TRWF 

2 


TRWR 

2 


fd^F^  ” Bd^F 


fd^R>  " Bdy*R 


(y0  - yu  ~ a * sin  fr) 

cos  ij; 

(y0  - yu  + b * sin  fr) 

cos 


for  yF  < 


for  y R < 


TRW 
2 ’ 

TRW 
2 ’ 


othe  rwise 


otherwise 


TRWF 

2 


TRWR 

2 


TRW 

2 

TRW 

2 


2.0  NUMERICAL  SOLUTION 

2.1  VERTICAL  TIRE  FORCES 

Given  zu,  4^,  9 ; z*u,  9 

Compute  axle  heights  ZT^,  ZT^  (1.5.2) 
Given  tire  force  vs.  deflection  fz  (ZT) 
Compute  tire  normal  loads  Fz^ 
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2.2  LATERAL  TIRE  SLIDING  FORCES 


Given  Fzi  , yu , r ; Compute  Fyi  (1.5.1) 

2.3  LATERAL  FORCE  AT  DISCONTINUITY 

Given  yu>  ip,  r ; Compute  yF,  yF>  yR,  yR  (1.5.3) 

• • 

Given  yF,  yF,  yR,  yR  and  discontinuity  force  vs. 
deflection  fd  (y.); 

Compute  FdF,  FdR , TRWF/2,  TRWF/2  (1.5.3) 

2.4  ROLL  CENTER  FORCES  AND  SPRUNG  MASS  ROLL  ACCELERATION 

• • • • 9 9 

Given  4>u  > ‘t’s » ^u>  ^s>  <ki>  ^s»  ^u’  ^s*  zs>  zu»  zs»  zu’ 

Compute  Spring  Deflections  and  Rates  (1.1.4) 

Given  Spring  Deflections  and  Rates,  fg 
Compute  suspension  forces  Fgg , FgR  (1.1.4) 

Given  <)>s , 4^,  ys,  yu,  zg,  zu,  4>s , Ys>  Yu»  zs>  zu>  FSL>  FSR> 
Compute  cj^  , Fl>  (1*1.1,  1.1.2,  1.1.3) 

2.5  COMPUTE  UNSPRUNG  MASS  MOTIONS 

Given  F^,  Fg^,  FgR,  4>s , tire  side  forces,  tire  normal  forces; 

Compute  yu,  zu,  ^ (1.2.1,  1.2.2,  1.2.3) 
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Given  tire  forces,  r,  4^; 

Compute  r,  0 (1.3,  1.4) 

2.6  UPDATE  VELOCITIES,  ORIENTATION,  AND  POSITION 

Given  0,  4>s  > 4\u  > ^ > Yu> 

Integrate  to  obtain  updated  states  and  rates 
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APPENDIX  B 


CREATING  THE  VEHICLE  CHARACTERISTIC  INPUT  FILE 


A.  OVERVIEW 

The  tripped  roll  over  simulation  program  is  designed  so  that  various 
vehicles  may  be  tested.  To  allow  the  user  easy  installation  of  each 
vehicle  to  be  tested,  a program  was  designed  to  create  an  input  file 
that  contains  the  necessary  physical  characteristics  of  each  vehicle. 
The  following  discussion  will  explain  how  J „ use  the  program  and  will 
define  each  of  the  characteristics  that  the  user  must  enter. 

To  run  the  program  from  DOS  the  user  types  either  TRFILECR  for  the 
unsprung  model  or  SPFILECR  for  the  sprung  mass  model  (underlined  phrases 
define  user  inputs)  while  defaulted  to  the  correct  disk  drive.  The 
program  will  initially  ask  for  a file  name.  The  name  should  include  the 
disk  drive  specifier  (if  necessary)  a file  extension  if  desired  and  a 
prefix  to  define  which  model  it  is  set  up  for  (e.g.,  A: STRNOVA. 75). 
Once  the  file  name  has  been  input  the  pro^  am  will  prompt  the  user  to 
enter  the  vehicle's  characteristics.  Examples  for  the  two  programs  are 
given  at  the  end  of  this  appendix.  Car  parameters  are  defined  in 
Table  B-l , and  values  for  a 1975  Chevrolet  Nova  are  given  in  Table  B-2. 

Before  any  parameter  file  is  created,  the  user  should  thoroughly 
read  Appendix  D — Parameter  Definitions  and  Specification.  Appendix  D 
gives  a descripiton  of  each  parameter,  including  the  units,  and  where 
they  can  be  found  in  the  equations  of  motion.  Furthermore,  Appendix  D 
gives  a detailed  explanation  on  how  to  find  some  of  the  uncommon  param- 
eters required  by  the  programs. 
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B.  TRFILECR  (Unsprung  Mass  Program  Parameters) 


ENTER  FILE  NAME  (TYPE  Q TO  QUIT) 


? TRNOVA.  7 5 


ENTER  THE  CAR  PARAMETERS  IN  THE  FOLLOWING  ORDER: 


ALPHA, BETA , TRW , HU , KZ , BZ , CF 
BD , KD1 ,KD2 ,MU , IX ,POINTl ,P0INT2 
IY , IZ , IY , WHLDEF ,KV ,TRAD ,HCURB 


SEE  PROGRAM  LISTING  OR  USER'S  GUIDE  FOR  REFERENCES  TO  THE  ABOVE 
TERMINOLOGY 


? 4.45,4.8,5,1.64, 12000,80, .75^ 

? . 25,24000, 96000, 121. 7,330, .16667, .6667  ? 
? 2240,2240,132,4000, 1,1, .5 g 

C.  SPFILECR  (sprung  mass  model  parameters) 

ENTER  FILE  NAME  (TYPE  Q TO  QUIT) 

? S TRNOVA  .15  2 


ENTER  THE  CAR  PARAMETERS  IN  THE  FOLLOWING  ORDER: 

ALPHA,  BETA,  TRW,  HU,  KZ , BZ,  CF , BD 

KD1  , KD2 , MU,  HR,  SUSK1,  SUSK2 , B1  , B2 

IX,  IY,  IZ,  HS,  IU,  WHLDEF,  KV,  POINT1 , POINT2 

MS,  IS,  TRAD,  HCURB , S,  SPRNGK,  SPRLNG,  ABSDEF , BL 


SEE  PROGRAM  LISTING  OR  USER'S  GUIDE  FOR  REFERENCES  TO 
THE  ABOVE  TERMINOLOGY 

? 4.45,4.8,5, .2,12000,80, .75,. 25  7 
? 24000,96000, 16. 7 , 1 , 6900,33000, 325, 33 2 
? 330,2240,2240,1.2, 132,4000, 1,. 1667, .6667^ 

? 105,240,1, .5, 1.75 ,200000, . 6667 , . 25, 1000  S 
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TABLE  B-l . CHASSIS  PARAMETERS,  CHEVY  NOVA  1975  MODEL 


ms 

= 

sprung  mass 

105  lb-sec^/ft 

hs 

= 

mg  c.g.  to  sprung  mass  bottom 

1.20  ft 

g 

= 

gravity 

32.2  ft/sec^ 

Xs 

= 

sprung  mass  roll  inertia 

240  lb-ft-sec" 

s 

= 

sprung  base/2 

1.75  ft 

K1 

= 

suspension  spring  constant 

6900  lb/ft 

B1 

= 

suspension  spring  damping 
constant 

325  lb-sec/ft 

B2 

- 

bump  stop  damping  constant 

33  lb-sec/ft 

K2 

= 

bump  stop  spring  constant 

33,000  lb/ft 

“u 

= 

unsprung  mass 

16.7  lb-sec^/ft 

TRW 

= 

track  width 

5 ft,  or  TRW/2  = 

= 

total  inertia  in  roll 

330  lb-ft-sec^ 

a 

= 

m c.g.  to  front  axle 

4.45  ft 

b 

= 

m c.g.  to  rear  axle 

4.80  ft 

Bz 

= 

total  inertia  in  yaw 

2240  lb-ft-sec^ 

= 

total  front  and  rear  unsprung 
mass  inertia  in  roll 

132  lb-ft-sec^ 

y 

= 

coefficient  of  friction 

0.75 

Kz 

= 

single  tire  rate 

12,000  lb/ft 

Bz 

= 

single  tire  damping 

80  lb-sec/ft 
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APPENDIX  C 


The  directory 
lowing : 


COMMAND 
SPTRIP87 
TRIP8  7 
SPTRIP 
TRIP 
SPFILECR 
TRFILECR 
SPFILECR 
TRFILECR 
GRAPHICS 
SPNOVA 
TRNOVA 
AUTOEXEC 
14 


Each  of  the  files 
familiar  with  all 


COMMAND.COM 


SPTRIP87.EXE 


TRIP8  7.EXE 


SPTRIP.  BAS 


TRIP. BAS 


DISK  CONTENTS 


(DIR)  of  the  tripped  roll  over  disk  contains  the  fol- 


COM 

17664 

3-08-83 

12: OOp 

EXE 

50464 

6-04-85 

10:36a 

EXE 

46336 

6-04-85 

2:  02p 

BAS 

35881 

6-04-85 

9:04a 

BAS 

'28988 

6-04-85 

1 : 36p 

EXE 

20040 

5-30-85 

8:57a 

EXE 

19624 

5-30-85 

8:49a 

BAS 

1114 

5-30-85 

8:53a 

BAS 

893 

5-30-85 

8:39a 

COM 

789 

3-08-85 

12 :00p 

75 

256 

3-26-85 

12:  04p 

75 

128 

3-26-85 

11:56a 

BAT 

78 

4-03-85 

7:56a 

File(s) 

106496  bytes 

free 

above  will  now  be  discussed  so  that  the  user  will  be 
the  files  on  the  disk. 


This  allows  the  computer  to  startup  and  load  the 
files  for  execution. 

This  is  the  compiled  version  of  the  program  that  is 
described  in  this  paper. 

This  is  the  compiled  version  of  less  complex  version 
of  SPTRIP.EXE. 

This  is  the  basic  interpreter  version  (code)  of  the 
program  described  in  this  paper. 

This  is  an  earlier  version  of  SPTRIP. BAS  without  the 
sprung  mass  degree-of-f reedom.  If  the  springs  in 
SPTRIP.  BAS  are  stiffened  so  that  the  sprung  and 
unsprung  masses  are  essentially  rigid,  the  results 
from  SPTRIP.  BAS  and  TRIP.  BAS  should  be  similar.  The 
reader  is  referred  to  Section  II-C,  Cases  21-26,  and 
Figs.  11-13.  Because  TRIP. BAS  is  less  complex  it 
will  run  faster  and  save  time. 
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SPFILECR.EXE 


This  is  the  compiled  version  of  SPFILECR.BAS  which 
will  be  discussed  below. 


TRFILECR.EXE 

SPFILECR.BAS 


TRFILECR. BAS 


GRAPHICS.COM 


SP NOVA.  7 5 


TRNOVA. 7 5 
AUTOEXEC.BAT 


This  is  the  compiled  version  of  TRFILECR. BAS  which 
will  be  discussed  below. 

This  program  is  used  with  the  SPTRIP  program.  It 
enables  the  user  to  create  a file  of  constant  car 
parameters  that  are  needed  in  the  SPTRIP  program. 
For  a more  detailed  look  at  this  program  the  reader 
is  referred  to  Appendix  B. 

This  program  is  used  with  the  TRIP  program.  It 
enables  the  user  to  create  a file  of  constant  car 
parameters  that  are  needed  in  the  TRIP  program.  This 
program  is  very  similar  to  SPFILECR.BAS  except  that 
less  inputs  are  needed  because  the  sprung  mass  is  no 
longer  present.  TRFILECR. BAS  can  be  used  by  follow- 
ing Appendix  B and  omitting  the  parameters  HR,  SUSK1, 
SUSK2 , Bl,  B2,  HS,  MS,  IS,  S,  SPRUNGK,  SPLNG,  ABSDEF , 
and  BL. 

This  program  allows  the  user  to  printout  the  time 

histories  that  the  SPTRIP  and  TRIP  programs  create. 
This  must  be  installed  before  SPTRIP  and  TRIP  are 

run. 

This  is  a file  of  the  constant  car  parameters  that 

were  used  in  testing  SPTRIP.  The  car  used  was  a 1975 
Chevy  Nova. 

This  is  a file  of  the  constant  car  parameters  that 

were  used  in  testing  TRIP. 

This  file  installs  the  graphics  made,  date,  time,  and 
heading  when  the  computer  is  booted,  or  switched  on. 
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APPENDIX  D 


PARAMETER  DEFINITIONS  AND  SPECIFICATIONS 

A.  INTRODUCTION  AND  OVERVIEW 

Based  on  experience  gained  by  using  the  "Tripped  Rollover"  programs, 
this  appendix  has  been  written  to  help  the  user  define  all  of  the  neces- 
sary parameters.  Table  D-l  lists  each  parameter,  their  units,  and  on 
what  page  in  the  Equations  of  Motion  Appendix  (Appendix  A)  the  parameter 
is  found.  In  addition  to  Table  D-l,  several  points  are  included  below 
that  will  provide  useful  guidance  to  the  user,  regarding  parameter  vari- 
ations. 

B.  MODIFIED  PARAMETER  DEFINITIONS  (TABLE  D-l,  FIGURE  D-l) 

1.  Sprung  Mass  Height 

HS  (hg)  plus  SPRLNG  add  up  to  the  sprung  mass  center  of  gravity 
height  above  the  axle  centerline.  The  value  for  SPRLNG  is  not  obtain- 
able from  static  tests  or  published  literature,  and  is  in  fact  a redun- 
dant parameter.  HS  can  be  set  to  equal  the  entire  distance  from  the 
c.g.  to  the  axle,  and  SPRLNG  is  then  set  equal  to  zero,  c.g.  height  is 
obtainable  from  static  test  data  or  published  literature.  A simple 
approximation  is  that  c.g.  height  = 0.39  x roof  height,  as  supported 

by  data  in  Fig.  D-2.  Referring  to  Fig.  D-l,  finally  note  that 

hs  ” ^c.g.  ~ ^R 
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2.  Equivalent  Spring  Rates 


SUSK1  (Kx),  SUSK2  (K2),  B1  (Bl),  B2  (B2)  should  be  set  to  their 

TRW 

equivalent  value  acting  at  the  wheels,  and  S (S)  should  be  set  = — j— , . 
The  reason  for  this  suggested  approach  is  that  the  main  effect  of  the 
springs  on  tip-over  is  their  effect  on  roll  stiffness.  The  roll  stiff- 
ness is  normally  obtained  from  static  tests  which  provide  a stiffness 
parameter  KrqlL  °f  moment/radian  of  roll  angle. 


In  a simple  case  of  a suspension  spring  = , this  roll  stiffness 

9 

would  be  equal  to  2K^S  . But  with  the  addition  of  auxiliary  roll  stiff- 

ness  devices  (anti-roll  bars),  the  roll  stiffness  is  higher  than  2K^S  . 

TRW 

When  setting  spring  rates  to  act  at  the  wheel,  S = -y  , then 

„ KX  (TRW)2 
2KlS2  = = KrOLL 

01  K = ^ KrQLL 

1 (TRW)2 


3.  Default  Values 

Several  parameters  in  the  model  are  primarily  included  to  minimize 
numerical  instability  problems  and  should  not  affect  rollover  potential. 
The  parameters  are  not  related  to  actual  physical  measurements,  and 
default  values  should  be  provided  as  follows: 

KV  = 1 

SPRNGK  = 200,000 
BL  = 1000 

C.  INITIAL  CONDITIONS 

In  order  to  provide  stable  initial  conditions  when  the  vehicle  is 
placed  in  a sliding  motion,  the  steady-state  roll  angle  and  tire  deflec- 
tions for  the  side-slide  conditions  should  be  entered  as  program  inputs. 
If  this  is  not  done,  the  vehicle  is  subjected  to  a step  input  lateral 
sliding  friction  force.  This  causes  the  vehicle  to  bounce  and  roll 
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which  in  some  cases  can  cause  the  vehicle  to  rollover  before  striking 
the  curb.  This  transient  may  also  lead  to  inaccurate  estimates  for  the 
minimum  lateral  sliding  speed  needed  for  tripped  rollover. 

The  equations  given  in  Table  D-2  should  be  used  to  calculate  initial 
conditions  for  <t>u,  <j)g,  Zu.  An  example  calculation  is  given  in 
Table  D-3.  Figures  D-3  - D-6  show  example  computer  runs  with  and  with- 
out the  corrected  values.  These  examples  were  run  in  the  initial  velo- 
city iteative  mode,  and  the  plots  shown  are  for  the  final  critical  side 
velocity  iteration. 

D.  NON-TRIPPED  TIP-OVEF  JNDITIONS 

A vehicle  can  tip-over  due  to  moments  generated  under  sideways  slid- 
ing conditions.  The  conditions  under  which  this  can  happen  are 
described  by  the  equation 

M2  Ms  gv  r Mu  g U r i 

Fzi  = f-  - c™r  [hs  + Tr1  ' 7™y  [Tr  + hJ 

where  F^,  -*-s  total  force  on  the  two  unloading  tires. 

Whenever  is  0,  the  vehicle  will  tip-over  during  sliding.  If 
the  program  is  run  under  these  conditions  in  the  velocity  iteration 
mode,  the  program  will  find  the  critical  side-slide  velocity  to  be  zero. 
By  setting  the  above  equation  to  zero,  the  c.g.  height  that  will  lead  to 
slide  induced  rollover  can  be  calculated.  Note  that  for  the  5 degrees- 
of-freedom  model  (i.e.  , jud  sprung  mass)  the  second  right  hand  side  term 
in  the  above  equation  must  be  dropped. 
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Figure  D-l . Revised  Parameter  Definitions 
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Center  of  Gravity  Height,  hc^(inches) 


35  - 
30  - 
25  - 
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PASSENGER  CARS 
O 1980  Chev.  Citation 
A 1975  Chev.  Nova 
0 Olds.  Cutlass  Ciera 
<$>  1956  Ford 
<]  Ford  Thunderbird 


• AMC  Eagle 
A Mitsubishi  Montero 
■ Suzuki  Samurai 
A Ford  Bronco  (full  size) 
^ Ford  Bronco  II 
^ Chevrolet  Blazer 


PICKUPS 
Nissan 
A Toyota 
Ef  Ford  F250XL 
<A  Chevrolet  SIO 

LARGE  VEHICLES 

(3  Dodge  Maxi-Van  School  Bus 
S Toyota  Dolphin  Motor  Home 
Ef  Chevrolet  Van  Ambulance 


Ar 


50 


60  70  80 

Roof  Height, h r00f  (inches  ) 


90 


100 


Figure  D-2 . Relationship  Between  Vehicle  Center  of  Gravity 
Height  and  Roof  Height 
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Figure  D-3 . 


1975  Neva  Test 


Case  at  Critical  Side  Velocity: 


Trimmed 
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Figure  D-4. 


1 Q7 5 Test  Case  at 


Critical  Ride  Velocitv: 


Not  Trimmed 
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Figure  D-5.  1975  Nova  Test  Case  Initial  Transients:  Trimmed 
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Figure  T-6.  1975  Nova  Test  Case  Initial  Transients:  Not  Trimmed 
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TABLE  D-l.  PARAMETER  DEFINITIONS 
FOR  TRIPPED  ROLLOVER  MODEL 


PROGRAM 

VARIABLE 

APPENDIX  A 
VARIABLE 

SEE  ON 
PAGE  NO. 

UNITS 

ALPHA 

a 

Horizontal  distance  from 
the  vehicle  center  of 
gravity  to  front  axle 

A-l 

ft 

BETA 

b 

Horizontal  distance  from 
the  vehicle  center  of 
gravity  to  the  rear  axle 

A-l 

ft 

TRW 

TRW 

Total  track  width 
(average  of  front  and 
rear) 

A-l 

ft 

HU 

hu 

Height  of  unsprung  mass 
center  of  gravity  above 
axle  centerline  (usually 
zero) 

A- 8 

ft 

KZ 

Kz 

Single  tire  radial  spring 
rate 

A-l  7 

lbs/f  t 

BZ 

Bz 

Single  tire  radial  damping 
rate 

A-l  7 

lbs -sec 
ft 

CF 

y 

Coefficient  of  friction 
between  the  tire  ground 

A-l  3 

— 

BD 

Bd 

Vehicle  structure  damping 
constant  in  the  lateral 
deformation  direction 

A-l  9 

1 bs  -s  e c 
ft 

KD1 

Wheel  and  suspension 
lateral  deflection 
spring  rate  when 
striking  the  curb 

A-l  8 

lbs/f  t | 

KD2 

Slope  of  force  vs. 
lateral  deformation 
of  body,  frame,  and 
engine  when  striking 
the  curb 

A-l  8 

lbs/ft 

MU 

Total  unsprung  mass 

A-l  0 

lbs-sec^ 

ft 
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TABLE  D-l.  (Continued) 


PROGRAM 

VARIABLE 

APPENDIX  A 
VARIABLE 

SEE  ON 
PAGE  NO. 

UNITS 

HR 

hr 

Height  of  the  lateral 
restraining  force 
member  (roll  axis) 
above  the  axle  center- 
line  (average  front 
and  rear) 

A- 8 

ft 

SUSK1 

K1 

Suspension  spring  rate 
per  side  equivalent 
acting  at  wheel 

A- 4 

1 bs  / f t 

SUSK2 

k2 

Bump  stop  spring  rate 
per  side,  equivalent 
acting  at  wheel 

A-4 

1 bs  / f t 

B1 

Suspension  shock 
absorber  damping  rate 
per  side,  equivalent 
acting  at  wheel 

A- 4 

lbs  -sec 
ft 

B2 

b2 

Bump  stop  damping  rate 
per  side,  equivalent 
acting  at  wheel 

A-4 

1 bs  -s  e c 
f t 

IX 

Xx 

Total  roll  inertia  about 
total  mass  center  of 
gravity 

A-l  2 

lb  ft/sec^ 

IY 

b 

Total  pitch  inertia 

A-l  2 

lb  ft/sec^ 

IZ 

Total  yaw  inertia 

A-12 

lb  ft/sec^ 

HS 

hs 

Height  of  sprung  mass 
center  of  gravity  above 
axle  centerline 

D-2 

f t 

IU 

Total  unsprung  mass 
roll  inertia  about  Mu 
center  of  gravity 

A-l  0 

lb  ft/sec^ 
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TABLE  D-l.  (Continued) 


PROGRAM 

VARIABLE 

APPENDIX  A 
VARIABLE 

SEE  ON 
PAGE  NO. 

UNITS 

WHLDEF 

— 

Constant  force  during 
crush  deformation  of  the 
"soft"  suspension  and 
wheel  members  when 
striking  a curb 

A- 18 

lbs 

KV 

*v 

Velocity  at  which  tire 
slides  with  a constant 
side  velocity  function 

A-14 

ft/sec 

P0INT1 

Deflection  distance  at 
which  wheel  and  suspen- 
sion begin  to  crush 
(permanent  deformation) 

A-18 

ft 

POINT2 

Crush  distance  at  which 
the  more  solid  members 
of  the  body,  frame,  and 
engine  become  involved 
in  lateral  deformation 

A- 18 

ft 

MS 

Ms 

Total  sprung  mass 

D-2 

lb  sec^ 
ft 

IS 

Sprung  mass  roll  inertia 
about  its  center  of 
gravity 

D-2 

lb  ft/sec^ 

TRAD 

tr 

Tire  radius 

D-2 

ft 

HCURB 

hcurb 

Curb  height 

A-10 

ft 

S 

s 

With  the  suspension 
spring  rates  set  to 
equivalent  acting  at 
the  wheels,  the  value 
for  s should  be  set  = 
TRW/ 2 

ft 

SPRNGK 

K 

Bushing  (lateral 
restraining  force) 
spring  constant 

A-8 

lbs/ft 
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TABLE  D-l.  (Concluded) 


PROGRAM 

VARIABLE 

APPENDIX  A 
VARIABLE 

SEE  ON 
PAGE  NO. 

UNITS 

SPRLNG 

— 

Redundant  term,  always 
set  = zero 

— 

— 

ABSDEF 

Vertical  wheel  travel 
distance  from  the  static 
position  of  the  wheel  to 
the  point  where  contact 
with  the  suspension  bump 
stop  begins 

ft 

BL 

B 

Bushing  (lateral  restraining 
force)  damping  constant 

A- 8 

lb-sec 

ft 
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TABLE  D-2 . INITIAL  CONDITIONS  TO  ACHIEVE  STEADY-STATE 
SLIDING  CONDITIONS  IN  THE  TRIPPED  ROLLOVER  PROGRAM 


<f\l  = 

MygM  [TR  + huJ  + Msgp  [hg  + TrJ 

Kz  (TRW)2 

^s  = 

MsgH  l-hs  " hrJ 

Kroll  + 

where  is  total  vehicle  roll  stiffness 

Initial  condition  for  static  tire  deflection 


Zu  - 

(Mu  + Ms  )g 
TR  + h»  ' 4 Kz 
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TABLE  D-3.  EXAMPLE  INITIAL  CONDITION  CALCULATIONS 


M = 16.7 
u 

g = 32.2 

y = 0.73 

Tr  - 1.0 

Ms  = 105 

h = 0.85 
s 

Kz  = 12,000 


TRW  = 5 


Xrqll  = 2Ki  S2  = 2 (3000)  (6.25) 


(16.7)  (32.2)  (0.75)  (1.0)  + 105  (32.2)  (0.75)  (1.85) 

12,000  (25) 

= 0.017  radians  (0.97°1 


<*>8 


105  (32.2)  (0.75)  (1.35) 

37,500  ^ 


0.108  radians  (6.2°) 


Z 


u 


1.0 


(121.7)  32.2 
48,000 


0.9185  feet 


Parameters  changed  from  Reference  test  case. 
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